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We develop a new approach to derive single state tight binding (SSTB) model for electron trans- 
port in the vicinity of valence-conduction bands of poly(G)-poly(C) and poly(A)-poly(T) DNA. The 
SSTB parameters are derived from first principles and are used to model charge transport through 
finite length DNA. We investigate the rigor of reducing the full DNA Hamiltonian to SSTB model 
to represent charge transport in the vicinity of valence-conduction band. While the transmission 
coefficient spectrum is preserved, its position shifts in energy. Thymine is poorly represented and 
its peak height is substantially reduced. This is attributed to the abstraction of the HOMO-LUMO 
coupling to other eigen-states in the nearest neighbor DNA bases, and can be corrected within 2""* 
order time independent perturbation theory. Inter-strand charge transport has also been analyzed 
and it is found that hopping to the nearest neighbor in the complementary strand is the most im- 
portant process except in the valence band of poly(G)-poly(C), where hopping to the second nearest 
neighbor between 3'— ends is the most dominant process. As a result, transport between 3'— ends 
and 5'— ends in the vicinity of valence band of polyG-polyC is asymmetric. 



I. INTRODUCTION 



II. METHODOLOGY 
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The growing interest in DNA as a molecular 
deviceii^ has led to a wide range of experimen- 
tal as well as theoretical work in the field. How- 
ever, charge transport results through DNA are still 

COntrOversial3i4i5,6,7,8,9,l"illd2,13d4il5,16,17,18,19,20,21^ 

perimentally it is found that DNA can either be 
a good conductor—, a semi-conductori^ or even an 
insulator—. Theoretically, while the exact charge trans- 
port mechanism is not clear, band transport 
polaronic transport ^^Ji fluctuation facilitated charge 
migrationjiSii^ and multi-step hoppingi^ii^ have been in- 
vestigated. The narrowing down of the precise charge 
transport mechanism has been difficult because the base 
pairs are weakly coupled due to the large inter-base sep- 
aration of ~ 3.4 A. This translates into relatively small 
hopping parameters for electrons between base pairs at 
the HOMO and LUMO energy levels (less than 125 mcV), 
which results in valence and conduction bandwidths of 
less than 500 meV. As a result of the narrow bandwidth, 
charge transport in DNA is easily modified by environ- 
ment effects, which include counter-ions, impurities, de- 
fects as well as hydration. Recently, ab initio calculations 
have focused on determining inter-base coupling^'^^'^*^ so 
that DNA can be represented within a single-state tight 
binding (SSTB) model. However, these derivations suf- 
fer from the limitations that they either cannot be ex- 
tended to describe interaction between different basest 
or can only find interactions between HOMO states (va- 
lence band) of DNA basesi^. Moreover, the rigor of re- 
ducing the full Hamiltonian to the SSTB model in de- 
scribing electronic transport through finite length DNA 
has not been addressed before. 

The aim of this work is to develop a general approach 
to obtain tight binding parameters to describe intra- 
strand and inter-strand interactions from first principles. 
This formalism will be applied to poly(G)-poly(C) and 
poly(A)-poly(T) DNA. 



We describe our model to determine the tight binding 
parameters for poly(G)-poly(C) here, and follow an iden- 
tical approach for poly(A)-poly(T). Initially one GO base 
pair is constructed using the Nucleic Acid Builder (NAB) 
software package22i2i24. The backbone of the structure 
generated by NAB is replaced by a hydrogen atoms, and 
the position of these H atoms are relaxed at the MP2 
level. The optimized structure of the hydrogen termi- 
nated DNA bases are shown in Fig^a. Using this con- 
figuration, double strand B-DNA structures'^ of four, six 
and eight GC base pairs are constructed. The self con- 
sistent Hamiltonian Hq of this structure is calculated 
using density functional theory, where the B3LYP den- 
sity functional and 6 — 31G basis set are used.^^ 
To work in an orthogonal basis set, we initially transform 
Ho to Hi such that. 



Hi = So-l/2HoSo-l/^ 



(1) 



where So is the overlap matrix. Every diagonal sub 
block of Hi, which corresponds to a DNA base, is diago- 
nalized and its eigenvectors are used to construct a block 
diagonal matrix U. Following this. Hi is transformed to. 



Ha = U^^HiU. 



(2) 



In this representation of H2, diagonal elements corre- 
spond to the localized energy levels of DNA bases and 
off diagonal blocks correspond to inter-base interactions. 
Independent from the simulated system size (four, six 
or eight base pairs), we find the hopping parameters be- 
tween HOMO and/or LUMO states beyond the second 
nearest neighbor base to be insignificant, and so they are 
neglected. Hence the Hamiltonian H2 is truncated and 
transformed to: 
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<"g,"g> 

ij = l^Afg 

+ $Z ^ns,j;"c j(C'ig,iCn^ J + C.C.) . (3) 

«Tig,,ie> 

£ng,i (crici) ^hc z*'' on-site energy of base ng (ric), where 
the subscripts g and c refer to guanine and cytosine re- 
spectively, tjig^i-n' ,j {tn^,i;n' ,j) IS the hopping parameter 
between energy levels i and j of base pairs Ug and n'g {ric 
and n^) respectively, tng^i-^n^j is the inter-strand hopping 
parameter between energy levels i and j of base pairs Ug 
and nc. iVg and Nc are the eigen-states in a single gua- 
nine and cytosine respectively. and C are the creation 
and annihilation operators and c.c. is the hermitian con- 
jugate. < ... > and ^ ... ^ represent 1** and 2"'' nearest 
neighbor interactions, respectively. 
In finding the parameters of the Hamiltonian in Eq. 
we have used the central two base pairs of the simulated 
system to minimize edge effects. The Hamiltonian in 
Eq. (PI is referred to as the full DNA model in the re- 
mainder of this paper. 

III. RESULTS 

A. SSTB parameters 

We present our results for on-site energy and intra- 
strand hopping parameters of the HOMO and LUMO 
states in table |l| In poly(G)-poly(C), the HOMO 
(LUMO) state is localized on guanine (cytosine), and the 
hopping parameter between consecutive guanines (cy- 
tosines) is 115 meV (61 meV). These values are much 
larger than the ones for poly(A)-poly(T), where the 
HOMO (LUMO) state is localized on adenine (thymine), 
and the corresponding hopping parameter between con- 
secutive adenines (thymines) is 21 nieV (23 meV). We 
have also determined the inter-strand hopping param- 
eters which are shown in Table ^1 Since the helical 
structure of the DNA breaks reflection symmetry in a 
plane perpendicular to the axis^^, we find that the inter- 
base hopping parameters depend on the directionality 
between 3'— and 5'— ends (Fig^b). The most strik- 
ing difference is for the HOMO state of poly(G)-poly(C), 
where tn <S: 3' - G - C - 3' >= 50 meV (dotted Une 
in FigHJb) and ii/ < 5' - G - C - 5' >= 7 meV (con- 
tinuous line in Fi^b). We find that the directionality 
dependence of the inter-strand hopping parameters cause 
significant asymmetry in inter-strand charge transport as 
discussed in the following sub-section. Finally, we note 



that both the intra-strand and inter-strand hopping pa- 
rameters do not depend on the system size as indicated 
by the results for the four and eight base pair systems 
shown in Tables HI and 



B. Transport Results 

Charge transport experiments typically involve either 
measuring the current-voltage characteristics of a DNA 
placed between metal contacts^ or measuring the charge 
transfer between donor and acceptor intercalators placed 
along DNA strands. While the parameters we derive 
here can be used to model both sets of experiments, we 
will focus on the former experimental configuration in the 
low bias limit. In presenting our results, we will compare 
the transmission probability in the vicinity of HOMO- 
LUMO states obtained using the full DNA (Eq. ©) and 
SSTB (Tables and models. However, to maintain 
the same injection rate from the leads to the device in 
comparing the two models, we account for all eigen-states 
of the edge base pairs that are coupled to the contacts. 
The transport calculations are carried out within the 
Landauer-Biittiker formalism^S and the transmission 
probability is, 

T = tr[r^GT^G°], (4) 

where C^"^ is the retarded (advanced) Green's function 
of the isolated DNA attached to the contacts, and 
pi-C^i) is the device coupling to the left (right) contact. 
While the coupling, F^^''^^ is crucial in determining 
the transport properties, it depends on the details of 
DNA-contact coupling, which is difficult to control 
experimentally. We have considered two limits of the 
DNA-contact coupling: the weak coupling limit where 
pL _ p7?, _ 20 meV and the strong coupling limit where 
pi _ pfi _ gQQ ]-neV^^. However, we present the results 
for the latter configuration. 

Finally, we note that in modeling intra-strand transport 
both DNA strands are coupled to the contacts at both 
ends. In the inter-strand representation only one strand 
is coupled to the contact at each end. 

Intra-strand transport: The results for intra-strand 
transport from the full DNA and SSTB models are shown 
by the solid and dashed lines respectively, in Fig. |2 
The SSTB model reproduces the peak and width of 
the transmission windows quite well, with the main dif- 
ference being a shift in both the HOMO and LUMO 
transmission windows to higher energies in the case of 
SSTB model. We find that this energy shift occurs 
due to coupling of the HOMO and LUMO states of a 
base pair to other eigen-states of neighboring base pairs. 
Further, because the hopping parameters between the 
HOMO /LUMO states of a base pair and the other eigen- 
states of neighboring base pairs is smaller than their en- 
ergy separation, we find that the energy shift seen in 
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Fig. |2] can be quite accurately accounted for using sec- 
ond order perturbation theory. The expressions for the 
second order correction to the HOMO (H) or LUMO (L) 
eigen-values are 

AE^'\ = - y Jk^B^. (5) 

We find that the energy shift for G (C) bases at the cen- 
ter to be —49 (—69) meV. The bases at the left and right 
edges have only one neighbor, and are shifted by —20 and 
—29 (—29 and —40) meV, respectively. When this shift 
is included, the eigenvalues from the SSTB model match 
the full DNA model. Calculations of the transmission 
probability which incorporate this second order correc- 
tion in the SSTB model agree more closely with the full 
DNA model as shown by the open circles in Fig. [21-a. 

The transmission probability for poly(A)-poly(T) are 
presented in Fig. \^h. It is clear from the right panel 
of Fig. El-b that the transmission window of thymine 
undergoes both a large energy shift and substantial 
peak height reduction compared to the case of poly(G)- 
poly(C). The reasons for this are: (i) weak intra-strand 
hopping parameters for the LUMO level of thymine 
compared to cytosine (table and (ii) large second 
order energy correction to the LUMO state due to strong 
interaction with other energy eigen-states. Following 
the perturbation theory analysis of Eq. (jsj, we find the 
second order energy correction to the LUMO eigen-state 
of thymine is —100 meV while that of the HOMO 
eigen-state of adenine is only —41 meV, for bases that 
are away from the edges. The corresponding corrections 
for bases at the left and right edges are —48 and —52 
(—16 and —25) meV respectively, for thymine (adenine). 
Again including these corrections in the Hamiltonian 
for the calculation of the transmission probability shows 
that the SSTB model for poly(A)-poly(T) agrees closely 
to the full DNA model as shown by the open circles of 
Fig.Hb. 

We now focus on the importance of including all 
eigen-states of the edge base pairs connected to the 
contact, in modeling charge transport. When the 
broadening due to the contacts and F-'^ is large, 
there is a non zero density of states in the valence 
and conduction band energy windows at the edge base 
pairs due to states other than en and e^. The hopping 
parameter between these states at the edges and, cr and 
tL of neighboring base pairs are non zero. As a result, 
charge injected into the tails of energy eigen-states 
other than and e^, contributes to the transmission 
probability at valence and conduction band energies. 
This contribution is referred to as the tail effect and can 
be partitioned out of the total transmission probability. 
In Fig.|21we present transmission spectrum between only 
HOMO (left panels of Fig. \^ and LUMO (right panels 
of Fig. [2Jl states at the contact as triangle symbols. 



Clearly, this corresponds to ~ 90% and 25% of the 
total transmission for G-C and A-T, respectively. It 
indicates that contribution of other eigen-states to the 
total transmission coefficient is important and they must 
be included by representing the contact with all modes. 

Inter-strand transport: Current-voltage measure- 
ments of inter-strand transport involve metal contacts 
connected to complementary strandsi^ That is, contacts 
are connected to cither only 5'— ends or 3'— ends of the 
DNA as shown in Fig. ^b. We first note that irrespec- 
tive of which strand an electron is injected into, transport 
occurs primarily along poly(G) (poly(C)) if electrons are 
injected energetically into the valence (conduction) band. 
Inter-strand hopping, which is the transmission limit- 
ing step, occurs mainly near the contacts because of the 
tiny density of states of the valence (conduction) band in 
poly(C) (poly(G)). We show the transmission probability 
for inter-strand transport in both the valence and con- 
duction bands of poly(G)-poly(C) in Fig. 01 The inter- 
strand transmission probability is more than an order 
of magnitude smaller than the intra-strand case because 
the density of states is peaked only along one strand in 
both the conduction and valence bands. Hopping into the 
strand with a smaller density of states limits the trans- 
mission / conductance in inter-strand transport. 

We will now gain some insight into the roles of the 
nearest and second nearest neighbor hopping parameters 
in table ITTl in determining inter-strand transport. This 
is done by calculating the transmission probability by 
setting specific inter-strand hopping parameters at the 
edges of the DNA to zero in the Hamiltonian. The solid 
(dashed) lines in Fig. |31-a,b correspond to setting the 
nearest (second nearest) inter-strand neighbor hopping 
parameter shown by the solid (dashed) line respectively 
in Fig.|31-c to zero. The solid triangles in Fig.|31-a,b is the 
reference, which corresponds to the full model. Clearly, 
in the conduction band, the nearest neighbor inter-strand 
hopping process shown by the solid line in Fig.|31-c is the 
most important in determining the transmission proba- 
bility between both the 3'— ends and 5'— ends. Further, 
the conduction band transmission probabilities across the 
3'— ends and 5'— ends are comparable. In contrast, in 
the valence band, the transmission probability across the 
3'— ends is twice as large as that across the 5'— ends. The 
reason for this is that the second nearest neighbor hop- 
ping parameter tn <^ 3' — G — C — 3' :s> determines the 
transmission probability across the 3'— ends, and is much 
larger (50 meV) than all other hopping parameters in the 
valence band of poly(G)-poly(C) as seen in Table ITU 
We have also investigated this asymmetric inter-strand 
transport for a hundred base pair (33.8 nm long) poly(G)- 
poly(C) system. We have found that it is is persistent and 
this represents a strong motivation to perform such inter- 
strand transport experiment. However, we note that 
these results are sensitive to the DNA base conforma- 
tion. Hence thermal effects can substantially reduce the 
asymmetric coupling along the DNA. Therefore, these ex- 
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periments should be performed at low temperature with 
ambient conditions that lead to B-DNA conformation. 
Inter-strand transport through poly(A)-poly(T) have 
also been carried out and showed symmetric transmission 
between the 3'— ends and 5'— ends. Inter-strand hopping 
occurs at nearest neighbor basis closest to the contact 
region. 

IV. CONCLUSION 



TABLE I: On-site energy and intra-strand hopping param- 
eters (meV) for HOMO (H) and LUMO (L) states. The pa- 
rameters are obtained from a four (eight) base pair system. 



Base 


£h 




+BB 


J.BB 


G 


-4278 (-4278) 


1137 (1148) 


19 (20) 


-115 (-114) 


C 


-6519 (-6533) 


-1065 (-1072) 


-61 (-60) 


-24 (-21) 


A 


-5245 (-5245) 


259 (258) 


24 (25) 


21 (21) 


T 


-6298 (-6346) 


-931 (-972) 


-23 (-23) 


-98 (-98) 



We have investigated the rigor of reducing the full 
Hamiltonian of a DNA to a single-state tight binding 
(SSTB) model. Tight binding parameters for both 
intra-strand and inter-strand transport have been tab- 
ulated. We find that the SSTB model quite accurately 
reproduces the transmission probability calculated from 
the full Hamiltonian, when second order corrections 
to the on-site potential of base pairs is included. As 
a result, the SSTB model is computationally efficient 
when compared to the full Hamiltonian. One caveat 
while applying the SSTB model is that injection of 
charge from the contacts to the edge base pair should 
be carefully modeled. This is because charge can be 
injected into the tails of all eigen-states of the edge base 
pairs, thereby making conduction due to the f,ail effect 
important. We have also investigated inter-strand charge 
transport and found strong asymmetric inter-strand 
current in the vicinity of the valence band (HOMO 
state) of poly(G)-poly(C). Current flowing between 
Ohmic contacts connecting the 3'-ends is almost twice as 
large as the current between Ohmic contacts connecting 
the 5'-ends. 



TABLE II: Inter-strand hopping parameters (meV) for 
HOMO (H) and LUMO (L) states. < ... > and < ... > cor- 
respond to nearest neighbor and 2"'' nearest neighbor inter- 
strand interactions, respectively. The parameters are ob- 
tained from a four (eight) base pair system. 



System 








tL 


tH 


<G-C> 








63 (63) 


2(5) 


<A-T> 








34 (34) 


26 (26) 


■C5'-G 


-C- 


5' 


> 


-12 (-12) 


-7 (-8) 


<3'-G 


-C - 


3' 


> 


-16 (-15) 


50 (48) 


« 5' - A 


~T~ 


5' 


> 


-10 (-10) 


-11 (-11) 


< 3' - A 


-T- 


3' 


> 


-13 (-13) 


9(9) 
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FIG. 1: (a) Atomic structure of hydrogen terminated DNA 
bases where arrows indicate H atom which replaces the back- 
bone, (b) 4 base pair G-C arranged in B-DNA configuration. 
5',3'-ends of DNA are shown and the sohd (dotted) arrow 
correspond to < 5' - G - C - 5' > (< 3' - G - C - 3' >) 
directional coupling. 



FIG. 2: Intra-strand conductance of (a) poly(G)-poIy(C) 
and (b) poly(A)-poly(T). Left (Right) panels correspond to 
transport through HOMO (LUMO) states. The linear re- 
sponse conductance is calculated at 300if using T{E) — 
J dE'T{E')f{E - E') where / is the Fermi function eval- 
uated at SOOif. Solid line - full dna model; Dashed lines - 
SSTB model with all eigen-states at the two edges; Open cir- 
cles - Same as "Dashed lines" but includes the second order 
energy correction given in Eq. 0; Solid triangles - full DNA 
model minus the tail effect. Insets: Same as main figures but 
conductance is calculated at bQK. 



FIG. 3: Inter-strand conductance (calculated at 300if) for 
poly(G)-poly(C). Left (Right) panels correspond to transport 
through HOMO (LUMO) states. In (a) only 5'- ends have 
an Ohmic contact while in (b) only 3'— ends have an Ohmic 
contact. Solid triangle - full model; Solid (Dashed) lines - 
Same as the solid triangle, except that the nearest (second 
nearest) neighbor shown by the solid (dashed) lines in (c) are 
set to zero. In (c), up-arrow (|) corresponds to the Ohmic 
contacts where charge is injected and collected. Each diagram 
in (c) corresponds to the graph in (a) and (b) directly above 
it. 
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